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PROCEEDINGS  OF  COLLOQUIUM  ON  STABLE 
SOLUTIONS  OF  SOME  ILL-POSED  PROBLEMS 


§1 .  Introduction 

Consider  the  following  operator  equation 

Af  =  g,  (0) 

where  A:  X  +  Y  is  a  mapping  from  a  Banach  space  X  into  a  Banach 
space  Y,  | f |  |g|  are  norms  in  X  and  Y  respectively.  If  the  inverse 
mapping  A  ^  exists  and  is  continuous  then  the  problem  of  finding 
feX  from  equation  (1)  is  called  well-posed.  In  many  cases  geR(A) , 
R(A)  is  the  range  of  A,  operator  A-1  exists  but  is  not  continuous 
and  is  defined  not  on  all  Y.  Then  the  problem  (1)  is  called  ill- 
posed,  because  f  does  not  depend  on  g  continuously.  In  what  follows 
we  assume  that  A  is  a  linear  bounded  operator  defined  on  X, 

Ker  A  »  {0},  geR(A) ,  A  ^  is  unbounded.  We  are  not  going  to  present 
here  the  most  general  results,  but  rather  concentrate  on  some  special 
situations  which  are  important  in  applications.  Another  feature 
of  this  presentation  consl sts in  the  following.  We  do  not  use  general 
methods  for  stable  calculation  of  the  solution  of  equation  (1)  such 
as  regularization  techniques  or  the  method  of  quasisolutions  ([1]— 
[5]),  because  these  methods  are  applicable  in  very  general  situations 
and  that  is  why  the  characteristic  features  of  some  special  classes 
of  problems  are  not  taken  into  account  in  the  process  of  applying 
these  methods.  On  the  other  hand,  these  general  methods  are  not 
easy  to  apply  from  the  standpoint  of  numerical  analysis .  The  methods 
used  in  what  follows  are  either  iterative  or  based  on  some  explicit 
analytical  formulas  for  approximate  solutions.  Yet  another  point 
to  be  emphasized  consists  in  the  following.  Suppose  that  the  mapping 
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A  is  compact,  X  *  Y  ■  H  where  H  is  a  Hilbert  space.  Operator  A~ 
is  unbounded.  Thus  the  problem  (1)  is  ill-posed.  But  there  are 
some  situations  when  it  is  natural  to  seek  the  solution  not  in  H 
but  in  a  wider  space  H_.  Operator  A  in  these  situations  can  be 
extended  to  the  operator  A:  H_  -*■  H+,  H+CHCH_,  the  embeddings  are 
dense  in  the  next  space,  and  the  extension  of  A  is  a  homomorphism 
from  H_  onto  H+.  Hence  we  made  of  an  ill-posed  problem  a  well- 
posed  one.  An  example  of  such  a  construction  is  given  in  §6.  The 
space  H_  in  §8  is  the  space  of  distributions.  Basic  problems  of 
the  static  fields  theory  can  be  reduced  to  integral  equations  of 
the  second  kind  at  a  characteristic  value.  This  is  another  kind 
of  ill-posed  problem  similar  to  ill-conditioned  problems  in  linear 
algebra  but  not  identical  to  them. 

In  §2  we  give  a  simple  remark  of  general  nature.  In  §3  an 
iterative  process  for  solution  of  integral  equations  with  positive 
kernels  is  given.  In  §4  an  iterative  process  for  solution  of  in¬ 
tegral  equation  of  interest  in  geophysics  is  given.  In  §5  an  iterative 
process  for  stable  solution  of  the  operator  equation 
f  =  X^  Af  +  g  will  be  given,  X^  is  a  characteristic  value  of  A. 

In  §6  a  study  of  integral  equations  basic  for  estimation  theory, 
filtering  and  prediction  of  random  fields  will  be  given.  It  is 
in  this  case  that  one  should  look  for  a  solution  in  the  space  of 
distributions.  In  §7  some  ill-posed  problems  are  discussed 
which  can  be  naturally  treated  as  the  problems  of  approximation 
theory.  In  §8  optimal  solutions  of  some  ill-posed  problems  will 
be  discussed.  In  §9  equation  with  many  solutions  is  treated. 

In  the  bibliographical  notes  we  discuss  some  other  trends 
in  the  theory  of  ill-posed  problems  and  give  references.  We  use 
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autonomous  numeration  in  each  section,  but  references  to  formulas 
of  other  section  have  the  number  of  section. 


This  introduction  we  end  with  examples. 

Example  1.  Let  F  be  a  perfectly  conducting  surface  in  not  neces¬ 
sarily  closed.  Then  equation 

/F  cr  (t)  (47re0rst)-:L  dt  =  1,  rgt  =  I  s— 1 1  (1) 

describes  the  equilibrium  charge  distribution  a  on  this  surface 
provided  that  the  potential  on  the  surface  is  equal  to  1,  eQ  >  0 
is  a  constant.  Equation  (1)  is  an  integral  equation  of  the  first 
kind. 

Example  2.  If  a  perfect  conductor  D  with  a  surface  r  is  placed  in 
the  electrostatic  field  -Vu  then  the  charge  distribution  a  on  T 
satisfies  integral  equation 

o  —  -Ac  -  2e0  3u/3N,  Act  =  Jp  cr(t)  2rrr —  (2) 

s  st 

N  is  the  outer  normal  to  r. 

It  is  known  that  -1  is  a  characteristic  value  of  the  operator 
A.  Also  it  is  known,  that  equation  (2)  is  solvable,  if 
/r  3u/3N  dt  =  0.  But  this  condition  is  satisfied  because  Au  *  0. 
Equation  (2)  is  so  called  an  equation  on  spectrum. 

Similar  equations  are  valid  in  other  static  fields  theories, 
e.g.  in  hydrodynamics,  magnetostatics,  elasticity. 

Example  3 .  Equation 

(x-t)-1  f(t)dt  *  g(x),  x  >  1  (3) 

is  of  interest  in  geophysics.  Its  multidimensional  analogue  is 
JD  f(t)  r”£  dt  »  g(x),  xeA,  AflD  »  0  (4) 

Example  4.  Basic  problems  of  random  fields  filtering  and  estimation 
theory  can  be  reduced  to  equation 
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(5) 


f.j  R(x,y)h(y)dy  =  f(x),  xeD  »  D  +  3u  ■  D  +  r, 
where  R(x,y)  is  a  correlation  function.  In  this  theory  the  solution 
of  equation  (5)  under  natural  assumptions  concerning  the  kernel 
R(x,y)  is  to  be  found  in  a  class  of  distributions. 

Example  5.  Some  problems  of  antenna  synthesis  and  apodization 
theory  can  be  reduced  to  equation 

fD  exp  (i(x,y)}  j(y)dy  ■  f(x),  xeA.  (6) 

In  optics  apodization  means  transformation  of  the  electro¬ 
magnetic  field  on  the  output  pupil  of  optical  instrument  in  order 
to  make  the  image  better,  to  improve  resolution  of  the  instrument. 
Example  6.  In  monoimpulse  radiolocation  the  direction-finding 
characteristic  (pelengation  characteristic)  is  defined  by  the 
formula 

g (k)  =  f2  (k+kQ)  -  f2(k-kQ),  (7) 

where 

f(k)«J*ff  j(x)  exp(ikx) dx  (8) 

is  the  radiation  pattern  of  a  linear  antenna,  j (x)  is  the  current 
distribution  and  0  <  kQ  <  it / ^  is  a  small  number.  Equation  (7) 
is  a  nonlinear  integral  equation  for  j(x). 

The  stable  approximation  to  the  solution  of  equation  ( 0) 
can  be  described  as  follows.  Assume  that  g  is  known  with  the 
error  <5  >  0,  i.e.  an  element  :  |g^  -  g|  5  is  known.  The  problem 

is  to  calculate  ffi  so  that  | f ^  —  f |  ♦  0  as  J  +  0  and  in  some  sense 
f ^  is  a  solution  of  the  equation  (o) .  For  example  one  can  require 
that  |Af5  -  g|  -*-0as5-»0. 

§2 .  Iterative  methods  of  stable  solutions  of  ill-posed 
problems . 


In  what  follows  we  omit  the  sign  of  space  in  the  notations  of  the 
norms  since  it  is  clear  what  norm  we  use.  Suppose  that  an  iterative 
process  is  known  for  solution  of  equation  (1)  ,  fR  =  Bn  g  -*■  f  =  A  ^g 
as  n  ■*  <*>,  and  the  operator  Bn  is  continuous  for  any  fixed  n.  Then 
it  is  possible  to  give  a  stable  approximation  to  the  solution  of 
equation  (1) . 

Indeed,  if  an  element  g^ ,  |g  -  g^ |  £  6  is  given  we  take 
^n(S)  =  Bn(6)  g<$  an<*  s*low  that  it  is  possible  to  choose  n(<5)  so 

that  If  -  f|  +  0  as  5  +  0.  To  this  end  we  start  with  the  in- 

n  ( 6 ; 


equality: 

I Bn96  "  flllBng  “  BngI  +  lBng  “  f I  •  (1) 

According  to  our  assumptions 

a(n)  =  I  Bng  ”  B  as  n  -*■  <*>,  (2) 

| Bng^  -  Bng|  <  b(5,n)  +  0  as  5  +  0  (3) 

for  any  fixed  n.  Thus 

lBng6  -  fl  1  a<n>  +  b(5,n) .  (4) 

Let  us  take  for  any  6  >  0  such  n  =  n(<S)  that 

a(n)  +  b(6,n)  =  min  =  a(S)  (5) 

From  (2) ,  (3)  it  follows  that 

n(j)  »  as  J  +  0,  a(j)  +  0  as  5  +  0.  (6) 

It  follows  from  (6)  and  (4)  that 

I Bn (5) g5  “  fl  -  0  as  5  -  0.  (7) 

We  proved  a  simple  but  useful  theorem. 


Theorem  1.  If  an  iterative  process  for  solution  of  equation  (1) 
is  known,  f  *  Bng  is  the  n-th  approximation  given  by  the  iterative 
process  and  Bn  is  a  continuous  operator,  then  it  is  possible  to  give 
a  stable  approximation  of  the  solution  of  equation  (1) . 
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n 


Remark  1.  If  B  =  7  ,  where  T  is  a  linear  bounded  operator, 

n  j-0 

then  |Bng5  -  Bng|  <  zMjT  5  3  b(5'n) 


Thus  in  this  case  there 


is  an  explicit  formula  for  b(6,n). 

Remark  2.  In  proposition  1  no  use  was  made  of  the  linearity  of  A. 


§  3 .  Iterative  process  for  the  solution  of  integral  equations  of 
the  first  kind  with  posihivg  kompls. 

Here  we  give  an  iterative  process  for  solution  of  integral 
equations  of  the  first  kind  with  positive  kernels.  Equation  (1.1) 
belongs  to  this  class  but  also  some  nonself-adjoint  equations  are 
in  this  class. 

Consider  the  equation 

Kf  =  / D  K(x,y)f(y)dy  =  g(x),  xeD  (1) 

r  2 

where  DCR  is  a  bounded  domain,  the  operator  K:  H  ♦  H,  H  =  L  (D) , 

is  compact,  K(x,y)  >  0.  Assume  that  a  function  h(x)  >  0  exists 

such  that 

Kh  <_  c,  J  a(x)dx  <  »,  a(x)  >_  m  >  0,  (2) 

where 

a (x)  s  h(x)  (Kh)"1  (3) 

Denote  by  <)>  =  f(x)a-1(x)  and  by  H+  the  spaces  L^(D,  a— 1(x)) 

2  2  xi 

with  the  norms  =  /D  jg|  a—  (x)dx.  Equation  (1)  can  be 

written  as 

K^  =  g,  K^  =  /D  K(x,y)a(y)<{>(y)dy.  (4) 

Let  K1  be  a  compact  operator  in  H+, 

Kl<f,j  *  Xj*j '  l  xil  >  U2l  !•••  (5) 
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The  kernel  in  (4)  is  positive.  We  suppose  that  some  itera¬ 


tion  of  is  continuous.  It  is  known  that  such  kernels  have  a 


unique  positive  eigenvalue  with  the  positive  corresponding 


eigenfunction  and  moduli  of  all  other  eigenvalues  are  strictly 


less  than  X^  (see  [12],  theorem  4.1.4).  We  rewrite  equation  (4) 


♦  *  Q*  +  g,  Q  =  I  -  Kr 


and  note  that 


Kx  Kh  =  Kh. 


Thus  operator  has  a  positive  eigenfunction  Kh  which  corresponds 
to  a  positive  eigenvalue  X^  =  1.  Hence  all  other  eigenvalues  of 
satisfy  the  inequality 


Aj|  <  1/  j  —  2 1  3  f  •  •  • 


Consider  the  iterative  process 


=  Q*  +  g,  =  g. 


Theorem  1.  Suppose  that  the  above  assumptions  about  K^(x,y)  hold, 
and  i)  equation  (4)  is  solvable  in  H+,  ii)  the  system  { 4> j }  forms 
a  Riesz  basis  in  H+  and  X^  ^  0,  |argX^|  <_  ir/3,  j  =  2,  3,.... 

Then  process  (9)  converges  in  H+  to  a  solution  <f>  of  equation 


(4),  equation  (1)  is  solvable  in  H_  and  its  solution  is  f  =  <f>  (x)  a-1  (x) 


Remark  1.  We  remind  that  a  system  { $ .}  forms  a  Riesz  basis  of  a 


Hilbert  space  H  if  the  system  is  complete  and  minimal  in  H  and  for 


any  numbers  C.,...,C  the  inequality 


a  I?lCj|2  <  lZ"c3i3l2  <  bIJ|Cj: 


holds,  where  0  <  a  <  b  are  constants  which  do  not  depend  on  n  (see 


[13]  )  . 
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Remark  2.  According  to  Proposition  1  iterative  process  (9)  can 
be  used  for  the  stable  approximation  of  the  solution  of  equation  (1) 


Remark  3.  In  applications  the  freedom  of  choice  of  a(x)  can  be 

useful.  For  example  the  solution  of  equation  (1.1)  for  an  unclosed 

surface  (thin  perfectly  conducting  screen)  must  have  the  known 

singularity  near  the  edge  L  of  F:  a~p  ^2(x),  where  p(x)  =  min|x-s| 

seL 

If  we  take  a(t)  =  p  ^'2  (t)  then  the  function  <j>  =  p^2(t)a(t) 
will  be  continuous  and  the  approximate  solution  of  equation  (1.1), 
obtained  by  means  of  iterative  process  given  in  Theorem  1,  will 
have  the  right  singularity  for  any  n.  It  is  clear  that  the  kernel 
of  equation  (1.1)  satisfies  conditions  of  Theorem  1. 

Remark  4.  In  [24]  iterative  process  (9)  was  applied  for  calculation 
of  the  magnetization  of  thin  magnetic  films.  This  problem  is  of 
interest  for  computer  technology. 

Proof  of  Theorem  1.  Let  6  =  $  -  <j>  .  Then  5  =  Qn6„.  Let 

n  n  n  u 

_  00 

<5q  =  This  expansion  is  possible  because  the  system 

{$.}  forms  a  Riesz  basis  of  H , .  We  have  6  =  TT  ( 1— X . ) nC . $ . , 
j  +  n  L1  j  j  j 

i5nlH*  1  b.  -  0  as  n  - 

Here  we  took  into  account  that  |argXj|  <_  ir/3  implies  |l-Xj|  <  1. 

Indeed,  if  X  =  rexp(i40,  r  <  1,  |^|  <_  it/3,  then  1 1-X  | 2  =  1+r2  - 
2 

2rcos4»  <  l+r  -  r  <  1.  We  proved  that  linuf  =  $  exists  in  H,  . 

—  n  + 

It  is  clear  that  <$>  is  the  solution  of  (4)  and  f  =  a<J>sH_  and  is 
the  solution  of  (1)  . 
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§4  Iterative  process  for  the  solution  of  an  itegral  equation  of 
interest  in  geophysics. 

Consider  the  equation 

(x-t)-1f(t)  =  g(x),  x  >  1.  (1) 

Its  multidimensional  analogue 

/Dl x-t | -1f (t) dt  =  g(x),  xe  A ,  AOD  =  0  (2) 

can  be  intepreted  as  the  problem  of  finding  mass  distribution  from 
the  known  potential. 

We  will  discuss  equation  (1)  in  detail  and  give  some  com¬ 
ments  concerning  equation  (2) .  Note  that  the  domain  of  integration 
and  the  domain  where  g(x)  is  known  are  different.  Let  us  set 
x  =  N+y,  -1  <_  y  <_  1,  g (N+y)  i  h(y)  ,  (3) 

where  M  >  3  will  be  chosen  later.  Equation  (1)  can  be  written  as 

Af  i  (N+y-t)-1  f  (t)dt  *  h(y)  ,  -1  <  y  <  1.  (4) 

This  equation  is  equivalent  to  equation  A*Af  =  A*h,  which  we  write 
as: 


Kf  =  t|»,  $  »  A*h,  K  =  A*A , 
where  the  kernel  of  K  is 


K  ( z ,  t) 


fl  _ dy 

J -1  (N+y-t) (N+y-z) 


1  on  (N+l-t) (N-l-z) 
t-z  (N+l-z ) (N-l-t) 


We  use  here  the  known  lemma  (see,  e.g.,  [7]). 


(5) 


(6) 


Lemma  1.  Let  A:  H  ■*  H  be  a  linear  operator  on  a  Hilbert  space, 
geD(A*)HR(A) .  Then  equation  Af=g  is  equivalent  to  equation 
A*Af  =  A*g. 


Proof  of  Lemma  1.  If  Af  =  g  and  geD(A*)  then  A*Af  =  A*g.  If 
A*Af  =  A*g  and  geR(A)  then  A*Af  =  A* Ah .  Thus  (A*A(f-h) , f-h)  =  0, 


2 

Operator  K  is  self-adjoint  and  nonnegative  in  H  =  L  [-1,  1] . 

Equation  (5)  is  equivalent  to  equation 

f  =  (I  -  K)f  +  ifr.  (7) 

From  (6)  it  follows  that 

K ( z , t) ~2N-2  as  N  +  (8) 

Hence  the  operator  B  =  I  -  K  is  positive  definite  for  large  N 
and  ||b||  <_  1.  Since  K  =  K*  is  compact  we  conclude  that  ||b||  =  1. 

Indeed 

1  |b 1 1  =  sup  |  f  -  Kf  [  >  |  <J>  -  K4>  j  =  1  -  X  +  1  as  n  +  ®,  (9) 

|f]<l  n  n  n 

Here  <j>n,  Xn  are  the  eigenfunctions  and  eigenvalues  of  K.  In  [  7] 
the  following  theorem  is  proved. 

Theorem  (M.  Krasnoselski j )  .  Let  B  =  B*,  S|b||  =  1  be  a  linear  operator 
on  a  Hilbert  space  H  and  -1  is  not  an  eigenvalue  of  B.  Then  the 
iterative  process 

fn+l  =  Bfn  +  g'  f0eH  (10) 

converges  in  H  to  a  solution  of  the  equation  f  =  Bf  +  g  if  this 

equation  has  a  solution. 

In  our  problem  operator  B  =  I  -  K  satisfies  the  conditions 
of  the  theorem.  Thus  we  have  proved  Theorem  1.  If  equation  (1) 
has  a  solution  in  H  then  iterative  process 

fn+l  =  +  f0£H  {11) 

converges  in  H  and  f  =  limfn  is  the  solution  of  equation  (1) 

provided  that  N  is  large  enough  (so  that  ||k||  <_  1)  . 

Remark  1.  Since  g  is  analytic  outside  the  cut  [-1,  1]  it  is 
uniquely  determined  everywhere  outside  this  cut.  The  function  f  is 
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uniquely  determined  from  the  jump  relations  for  the  Cauchy  integral 

Remark  2.  According  to  Theorem  1  from  §2  iterative  process  (11) 
gives  a  possibility  to  construct  a  stable  approximation  to  the 
solution  of  equation  (1) . 

Remark  3.  The  method  of  solution  of  equation  (1)  holds  for  equation 
(2)  without  essential  alternations,  but  the  question  about  uniqueness 
of  solution  of  equation  (2)  is  more  difficult.  Under  some  addi¬ 
tional  assumptions  about  the  solution  of  equation  (2)  uniqueness 
of  the  solution  can  be  proved  [14]. 

From  the  Krasnoselski j  theorem  the  following  proposition 
follows.  Let  Af  =  g,  A  is  a  bounded  linear  operator  on  H,  geR(A) . 
Then  the  iterative  process 

fn+l  "  ~  aA*A)fn  +  aA*g,  fQeH,  0  <  a  <  [J^> A 1 1  <12) 

converges  in  H  to  a  solution  of  equation  Af  =  g  for  any  fQeH. 

§5.  Iterative  process  for  the  solution  of  the  equation  of  the 
second  kind  at  a  characteristic  value. 

If  a  perfect  conductor  is  placed  in  the  initial  field  with 
the  potential  -Vu  then  the  charge  distribution  on  the  surface  of 
the  conductor  satisfies  equation  (.1.2).  The  number  -1  is  the 
smallest  characteristic  value  of  the  operator  A  in  equation  (1.2) 
and  this  equation  is  solvable  since  /p  3u/3N  dt  =  0.  Since  equation 
(1.2)  is  an  equation  at  a  characteristic  value  (equation  on 
spectrum  we  shall  write  in  the  sequel) ,  the  problem  (1.2)  is  ill- 
posed.  Similar  equations  occur  in  hydrodynamics  of  ideal  incom¬ 
pressible  fluid  and  in  elasticity.  In  all  these  cases  the 
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characteristic  value  is  semi-simple.  It  means  that  the  root  space 
corresponding  to  this  characteristic  value  coincides  with  its 
eigenspace.  We  present  now  the  abstract  setting.  Let  A  be  a 
linear  operator  on  a  Hilbert  space  H.  Suppose  that  the  spectrum 
of  A  is  discrete,  Aj  are  the  characteristic  values  of 
A,  <|>j  =  A  j A«f>  j ,  |  AjJ  <  j  A ^  I  £  |  X ^  I  <_....  Suppose  that  A^  is  semi¬ 

simple.  Denote  by  G  =  Ker  (I  -  T^A*) ,  by  { ^ j > j = i  an  orthonormal 

basis  of  G,  by  G"^~  =  R(I  -  A^A),  where  R(A)  is  the  range  of  A,  by 
P  orthoprojector  in  H  onto  G.  We  consider  the  equation 

g  =  A.Ag  +  f,  feG^-,  ^  (1) 

S’ 

and  introduce  the  following  operator 
m 

Bvg  =  Ag  +  y  l  (g,t|>.)if» ..  (2) 

Y  j=l  :  D 


Denote  by  r 

Y 

r  =  min  ( ! A 

Y  1 


the  following  number 
2\,  |A1(1  +  Y^)’1!  )  , 


(3) 


where  y  is  a  numoer  which  we  fix  later. 

Equation 

9  *  xiV +  f  j_  <4> 

is  equivalent  to  equation  (1)  on  the  set  G  .  It  means  that  every 
solution  geG  of  (1)  is  a  solution  of  (4)  and  vice  versa.  Let 
us  assume  that  dim  G  •  dim  ker  (I  -  A^A) .  Tnis  is  so  if  e.g.  A  is 
compact. 


Theorem  1.  Operator  B,  has  no  characteristic  values  in  the  disk 


| X |  <  r  .  If  i 1  +  yX^\  <  1  then  the  iterative  process 
9n+l  "  xlVn  +  F'  90  *  F' 


(5) 


converges  not  slower  than  geometrical  progression  with  denominator 
q,  0  <  q  <  |AjJ  r  *  to  the  solution  of  the  equation  g  =  A^B  g  +  F, 


VFeH. 
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If  FeG  then  F  can  be  represented  in  the  form  F  =  X^Av  -  v  where 
v  is  some  element  of  H.  In  this  case  g  =  <J>  -  v,  with 

Ker  (I  -  X^A)  and  P4>  =  Pv.  If  dim  G  =  1,  <f>e  Ker  (I  -  X^A)  , 

||<Ji  ||  =  ||<M|  =  1'  then  $  =  $(£,$)/($,$).  Process  (5)  is  stable  in 
the  following  sense:  the  sequence 

Vl  ■  XlByhn  -  F  +  En-  V  F'  M  i  E  l6) 

satisfies  the  estimate 

lim  sup  llg  "  h  ||  «  0(e)  as  n  (7) 

where  g  =  lim  g_  as  n  ■*  ». 
n 

Theorem  2 •  The  iterative  process 

£n+l  -  XlA£n'  £0  *  £  181 

converges  at  the  rate  of  geometrical  progression  with  denominator 

q,  0  <  q  <  | X1X~1 |  to  an  element  of  Ker  (I  -  X^A) .  If  dim  G  =  1 

then  this  element  is  equal  to  <j»(f  ,<p)  /  (<p  ,ip)  . 

Remark  1.  The  problems  of  static  fields  theory  (electrostatic, 
hydrodynamics,  elasticity)  lead  to  equation  (1) .  In  electrostatics 
(equation  (1.2))  dim  G  =  1,  in  elasticity  and  hydrodynamics  dim  G  = 
=6.  In  electrostatics  when  the  boundary  r  has  n  connected 
components  (the  case  of  n  bodies  in  the  exterior  field)  dim  G  =  n. 

Proof  of  Theorem  1.  Let  g  =  XB^g.  Multiplying  this  equation  by 

we  get  (g,  <(>..)  =  XX'^g,  i|k)  +  Xy(g,  <J^j).  Thus 
✓ 

(g,  ’J'j)  (1  -  xx^1  -  Xyj  =  0.  If  (g,  t|ij)  ?  0  for  some  j,  1  <_  j  £  m, 

then  X  =  X^  (1+yX^)-1.  If  (g,  ip.)  =  0  for  1  £  j  <_  m,  then 

B^g  *  Ag  and,  therefore,  Xec(A)  (o (A)  is  the  spectrum  of  A).  From 

this  we  conclude  that  in  the  disk  I X I  <  r  ,  where  r  is  defined 

Y  Y 

by  formula  (3),  there  are  no  characteristic  values  of  Indeed, 
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if  X^  ea(B^)  then  (g,  ^  j )  ■  0  for  g^O  and  hence 

g  *  X^Ag,  geG  .  This  is  impossible  because  X^  is  semisimple. 

2 

(We  remind  the  reader  that  X^  is  semisimple  iff  (I  -  X^A)  g  =  0  => 
=>(I  -  X^A)  g  =  0,  or  iff  {(I  -  X^A)  g  -  0,  geG  }  =>g  »  0,  or 
iff  X ^  is  a  simple  pole  of  (I  -  XA)”1).  Hence  X^fafB^).  Spectrum 
of  B^  consists  of  ofAJXX^^  and  possibly  of  number  X^l  +  yX^)”  . 

If  $  =  XB^,  X  ^  X^,  then  4>eG  ,  J>  =  XAt  and  therefore  Xec (A) . 

Let  us  prove  that  process  (5)  converges.  If  we  choose  y 


-1 


so  that  |1  +  yX^|  <  1,  then  q  =  |X^!r^  <  1 

there  are  no  characteristic  values  of  B  in  the  disk 


and  we  conclude  that 


<  r. 


Thus  process  (5)  converges  not  slower  than  geometrical  progression 

-1 


with  denominator  q  = 

1 


X]_|r 


0  <  q  <  1.  If  F  *  A.  Av  -  v  it  is 
00 


clear  that  FeG  .  Thus  Af  =  B  F.  Let  g 


I 

i 


l  X?  B3F,  then 

j-0  1  Y 


geG  ,  Ag  =  B^g,  g  *  X^B^g  +  F.  From  this  it  follows  that 

g  +  v  *  XjA(g  +  v) ,  so  that  g  +  v  *  te  Ker  (I  -  X^A) .  Since 

Pg  =  0  we  conclude  that  Pv  =  P$ .  if  dim  G  =  dim  Ker  (I  -  X^A)  =  1, 


(j>eG,  ||4  [|  =  i|tj|  =  1,  then  $  =  c$,  c  =  const.  Multiplying  this 
equality  by  if)  we  obtain  c  *  ($,t)/(f ,ip) ,  because  (g,ii>)  =  0  and 
(<t>,<p)  ?  0.  It  remains  to  prove  the  stability  of  process  (6), 
i.e.  estimate  (7) .  We  have 

h„-  jo  (xiV3  F  +  X  ‘W1 


9  -  l  ‘W3  F'  llXlByll  i  9' 
j=0  '  ' 

From  this  it  follows  that 


n-1 


;g-hnll  £  l  <I3e  +  l  qJ||F||<  (e  +  |jF  ||qn+1)  (1-q) _1 . 
j-0  j*n+l 


(9) 
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Formula  (8)  follows  from  (9) 


Proof  of  Theorem  2.  First  we  prove  the  following  lemma. 

Lemma  1.  Let  a  function  F(z)  be  analytic  in  the  disk  |zf  <  r 

and  meromorphic  in  the  disk  |z|  <  r  +  e,  e  >  0  with  values  in  the 

algebra  of  linear  bounded  operators  on  a  Hilbert  space.  Suppose 

that  X  is  a  simple  pole  of  F(z),  res  F(z)  *  c,  |X|  =  r.  If  there 

z=X 

are  no  other  poles  of  Fiz)  in  tne  disk  [ A |  <  r  +  z  and 

F(z)  =  T  a  zn  for  1 z ■  <  r,  then  -  lim  Xn+  a  -  c. 

n  11  n 

h=0 

Proof  of  Lemma  1.  The  function  F(z)  -  c(z  -  X)  *  is  analytic  in 

the  disk  |z|  <  r  +  e.  Hence  F(z)-c(z-X)~  =  I  b  zn, 

n=0  n 

oo 

|zi  <  r  +  e.  For  jzj  <  r  the  identity  holds:  £  b  zn  = 

n»0  n 

*  7  (a  +  cX  i)zn.  This  identity  is  valid  in  the  disk 

n=0  n 

|z|  <  r  +  e  according  to  analytic  continuation  principle.  Thus 

lim  (Xn+1  a„  +  c)  *  0. 
n~  n 

Passing  over  to  the  proof  of  Theorem  2  we  note,  that  the 
,  «  .  . 

function  (I  -  zA)  ”  f  »  £  z-*A-*f  is  analytic  in  the  disk 


{ z l  <  X^,  has  a  simple  pole  at  z  *  X^  and  has  no  other  poles  in 

the  disk  |zj  <  X2.  Applying  Lemma  1  we  conclude  that 

lim  (x”+^  An+^  f  +  c)  *  0.  The  rate  of  convergence  is  0  (|X^  X^i”) 

Since  f  ■  X?  Anf  we  have  lim  f  =  g  exists.  It  is  clear  that 
n  1  n 

n-*-00 

g  »  X^  Ag.  If  dim  Ker  {I  -  X^A)  *  1  we  have  g  *  a$,  a  =  const  and 
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m 


♦  e  Ker  (I  -  X.A),  ♦  +  0.  We  note  that  4>)  =  (X.  Af  ,  t|i)  « 

1  n+l  1  n 

■  (£n,  4>)  (f,  <J>)  •  Thus  a  (<(>,<]>)  =  (f,<|>),  a  =  (f  ,$)/($  ,ip)  , 

Remark  2.  Theorem  2  gives  an  abstract  analogue  of  the  iterative 
process  for  finding  the  equilibrium  charge  distribution  on  the 
surface  of  a  perfect  conductor. 

§ 6 .  How  to  make  of  an  ill-posed  problem  a  well  posed  one. 

In  general  there  are  several  answers .  One  way  is  to  substitute 
the  original  problem  Af  *  g,  A:  H  ■*  H  by  some  other  problems 
(methods  of  quasi  solutions  and  regularization) .  The  other  way 
is  to  consider  an  extension  of  the  operator  A,  A:  H_  ■»  H+, 

H+  c  H  C  H_  and  to  choose  H_  and  H+  so  that  A  will  be  a  homomorphism 
of  H_  onto  H+.  In  this  case  equation  Af  =  g  has  a  unique  and 
stable  solution  in  H_.,  f  ®  A  *g.  We  give  an  example  of  a  class 
of  integral  equations  which  can  be  treated  in  such  a  way.  The 
basic  integral  equation  in  random  fields  filtering  and  estimation 
theory  can  be  written  in  the  form 

fD  R(x,y)  h(y)dy  =  f(x),  x£DCRr,  (1) 

where  D  is  a  bounded  domain  with  a  smooth  closed  boundary  ", 

D  »  D  +  r,  r  >  1. 

Following  [16]  let  us  assume  that  the  kernel  R(x,y)eR. 

2  r  2  r 

It  means  that  a  self  ad  joint  elliptic  operator  L:  L  (  R  )  L  (  R  ) 

and  polynomials  P(X)  >  0,  Q(X)  >  0,  Xel  *  (-«,<*>)  exist  such  that 

R(x,y)  -  /A  P (X ) Q-1 ( X ) $ (x,y , X) dp (X)  (2) 

Here  A,  <5,  dp  are  the  spectrum,  spectral  kernel  and  spectral 
measure  of  L  respectively.  Let  p,  q,  s  be  deg  P,  deg  Q  and  ord  L 
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respectively,  and  a  =  0.5(q-p)s.  Denote  by  the  scale  of  Sobolev 

spaces  (D) ,  H_fc  for  t  >  0  is  the  dual  space  to  Hfc  with 

2 

respect  to  Hg  *  L  (D) .  It  was  proved  in  [16]— [19]  thatthe  operator 
R:  H  +  H  is  a  homomorphism  of  H_a  onto  H^.  It  means  that  the 

solution  of  equation  (1)  in  H_a  exists,  is  unique  and  stable 
towards  small  perturbations  of  f  in  Ha. 

Some  explicit  analytical  formulas  for  the  solution  of  equa¬ 
tion  (1)  with  kernel  R(x,y)eR  were  obtained  in  [17].  This  leads 
to  the  following  method  of  stable  solution  of  some  integral  equations. 
Given  a  kernel 

R(x,y)  =  /A  R(X)*(x,y,X)dp(X) ,  (3) 

where  R(X)  is  a  function  which  satisfies  the  following  condition: 

0  <  R(X)~A(1  +  X2)"*8,  S  >  0  is  integer.  (4) 

Then  for  any  e  >  0  there  exist  polynomials  P£(X),  Q£(X),  deg  Q£(X)  - 
-  deg  P£(X)  =  26  such  that 

sup  {(1  +  X2)8  |R(X)  -  R  (X)  | }  =  ||R(X)  -  R  (X;  ||  <  e,  (5) 

Xe  I 
where 

R_ (X )  =  P  (X)  Q”1 (X ) ,  (6) 

c,  Z  Z 

and  I  =  (-00,00)  .  Consider  the  equations 

Rh  -  f,  R£h£  =  f,  (7) 

where  the  kernels  of  the  operators  R  and  R£  are  defined  by  the 
functions  R(X)  and  R£(X)  respectively. 

Theorem  1.  If  (4 )  holds  then  for  any  there  exist  unique 

solutions  h,  he  H  -  of  equations  (7)  and 

£  "(JS 

I  h  -  h£  |  _gg  <  2M2e(l-2eM)_1|f|eg,  (8) 
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provided  that  2eM  <  1,  where 

M  -  rinf  inf  [(1  +  X2) S | R  (X ) | ] ,“1.  (9) 

le>0  Xel 

Proof  of  Theorem  1.  We  start  with  the  following  lemma. 


Lemma  1.  Let  H+  c  H0  C  H_  be  a  triple  of  Hilbert  spaces,  H_  is 
the  dual  space  for  H+  with  respect  to  HQ.  Suppose  that  R:  H_  -*•  H+ 
is  a  linear  map  and 

cx!h;f  <  (Rh,  h)  <  cl  hi  l,  cx  >  0, VheH_  •  (10) 


Then  R  is  a  homomorphism  of  H_  onto  H+, 

i  R,  <_  2c,  |  R  i  £  c^  ,  ( 

where  ,R;  is  the  norm  R:  H_  -*•  H+,  jhi_  is  the  norm  in  H_. 

Proof  of  Lemma  1.  Operator  R  is  defined  on  all  H_  and  from  (10) 
it  follows  that 


iRh,+  >  c^h^ 

Therefore  R*1  is  defined  on  imR  =  range  R  and 


(12) 


(13) 


It  remains  to  prove  the  first  inequality  (11) .  We  have 
j  R j  *  sup  |  (Rg,  h)  |  *  sup|  ^  {  (R(h+g)  ,  h+g) 

i»:.H 

-  (R(h-g)  ,h-g)  -  i[  (R(h+ig)  ,h+ig)  -  (R(h-ig)  ,h-ig)  ]  }  |  <_ 

<_  sup  |  +  |h-g|2  +  |h+ig|2  +  jh-ig|2} 

i  g i  _<l » I h | _<_  1 

<  c  sup  { j h j 2  +  j gj 2}  <  2c  (14) 

Pi.11 

It  remains  to  prove  that  imR  ■  H+.  This  is  true  because  R  is 
monotone,  coercive  and  continuous  (see  [20]). 
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Lemma  2.  Suppose  that  R^ :  H_  -*•  H+  is  a  linear  bijective  mapping 
of  H_  onto  H+  for  any  e  >  0,  0  <  e  <  Eg,  where  >  0  is  some 

fixed  number.  Suppose  that  R:  H_  4  H+  is  a  linear  mapping, 

^Re^|  £  m ,  |R  -  R  |  <  e,  where  M  >  0  does  not  depend  on  e, 

0  <  e  <  Eg.  Then  R  is  a  bijection  of  H_  onto  H+  ,  and 

|R-1|  <  M  (1  -  eM) ~1  for  eM  <  1  (1! 

Proof  of  Lemma  2.  We  have  R  =  R£  [  I  +  R^  (R  -  Rp)],  where  I  is 
the  identity  in  H_.  The  operator  R*1  (R  -  R£)  acts  in  H_  and 

Ir"1  (r-r  )  |  <  cm.  If  cm  <  1  then  R-1  ■  [I  +  r'1  (r  -  R=.)]“1r"1 


Ir"1!  <  m  (i  -  eM)"1 


Now  it  is  possible  to  prove  theorem  1.  Using  Parseval  equality 


we  conclude  that 


(Rh,h)  *  JA  R(X) |h(X) l^dq(X)  <  c  JA  (l+X2)~Blh!Zdq(X)  =  c|h|jgs.  (17) 

Here  h  is  the  image  of  h  in  the  representation  generated  by  the 
expansion  in  eigenfunctions  of  the  operator  L,  and 

c  =  sup  { (1+X2) 3  R (X ) | } .  (18) 

Xel 

If  we  denote 

c.  -  inf  { (1+X2) 3  R( X ) ) ,  (19) 

1  Xei 

then 

C1  )h)fSs  £  (Rh'  h>  £  cihLgs-  (20) 

Prom  (20)  and  Lemma  1  it  follows  that  R  is  a  bijection  of  H_gs  onto 
Hgs  and  (11)  hold. 

From  (9)  and  Lemma  1  it  follows  that 

I  Re1!  £  M.  (21) 
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From  (5) ,  (18)  and  (11)  it  follows  that 

|R  -  Rel  <  2e  (22 

From  (21) ,  (22)  and  Lemma  2  it  follows  that 

|  R”1 1  <_  M  ( 1  -  2eM) _1  for  2eM  <  1.  (23 

We  have 

|R-1-R~1|  =  Ir"1  (R£-R)R-1|  <  | r”1 |  | R-Re 1  i R-1 1  1  2£M2(l-2eM)'1. 

Inequality  (8)  follows  from  (24)  .  This  completes  the  proof 
of  Theorem  1. 

§7.  Some  problems  of  approximation  theory. 


1.  Suppose  that  f(x)eC^  is  a  given  function,  =  C^(I), 

I  =  [-1,  1],  ||f|j  =  max  |  f  (x)  |  ,  E  (f)  =  min  ||f(x)  -  P(x)||, 

xel  P (x) eP 

n 

?n  is  the  set  of  all  polynomials  of  degree  <_  n,  w(f;t)  = 


=  sup 
I  x-y 


f(x)  -  f  (y)  |  ,  Tn(x)  =  cos  (n  arc  cosx)  , 


<^t,x,yel 


Xj  *  xj , n+1  *  cos  fjvFI  1  <  j  <  n+i, 


T  , . (x) 

1  (x)  «  l.(x)  -  - — - 

3  3  TA+l(xj) (x_xj} 


h+1 


L  (x;  f )  =  I  f<*3>VX)' 
n  j-1 


(1) 


LnS>(x'fi»  * 

n+1 

jii  f«,x5> 

n+1 

Xn  " 

n.i,  hj 

3=1 

(x)  |  ||  , 

h+1 

Xn(x) 

-  1  1 

j-1 

lj(x) 1 , 

X’  » 
n 

l|x;(x)  !|. 

(2) 

(3) 

(4) 
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(5) 


r 


mpmmm 


lunji.pii  m 


Given  a  set  of  numbers  f^(Xj),  1  £  j  ±  n+1,  jf^(Xj)  -  f(x^)|  £  5 

we  are  going  to  find  a  stable  approximation  for  f(x)  and  give  an 
error  estimate  of  the  formula  for  approximation. 


Theorem  1.  If  feC^  then 

||f  (x)  -  Ln(x,f)  ||  <  (l+\n)En(f)  , 

I |f  *  ( x)  -  I^(x,f)||<  (l+^n)En-i(f ' 

nX 

if'(x)  -  L'  (x,  f )  |  <  (1  +  - — )  E  ,  Cf '  )  , 

n  t  ii""  x 

/1-x2 


u  —  ri  11  —  - 

/1-x2 

||f  |x)  -  l/k)(x,f)||  <_  C  wfn"1;  f(r))nk_r  (l+nklnn)  ,  k  <_  r. 


(6) 

(7) 

(8) 


(9) 

(10) 


Theorem  1'.  If  feCr,  r  >  3,  then  there  exists  n(5)  such  that 


llLnm  (X'V  “  =  6  ^  °'  (11) 

[]dLn(|l!X>V  -  f •  ||  =  0  (Sr~2  IMS-1)),  5-0  (12) 

i|x  Ln(J)(x'f5)  -f’l  iEn-l(f’)(1+XA(x))  +5XA(x)'  (13) 

(l+Xn)En(f)  +  SXn,  (14) 


Proofs  of  theorems  1  and  2  are  given  in  [21].  The  following 
Lemma  is  used  in  the  proof  and  also  is  of  independent  interest. 
Lemma  1.  Suppose  that  P^(x)£Pn  and 
m 

I  |P, (x) |  <  M,  xel,  (15) 

j-1  2 

where  m  is  arbitrary.  Then 
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m  imp  vpiipni  wur 


m  Mn  2 

l  |pl(x)|  <  min  (  . - *•  ,  Mn  )  .  (16) 

j=l  J  ~  /l-x 


2.  The  following  approximation  problem  is  of  interest  in  optics 

and  electrodynamics  ([19], [29]). 

Let  DCRr  and  ACRr  be  some  simply  connected  and  connected 

2 

domains,  f(x)eC(A)  or  f(x)eL  (A).  Denote  by  the  class  of  func¬ 
tions  which  can  be  represented  in  the  form 

g(x)  *  (2tt)”N//2  Jd  exp(-i(x,y)  )h(y)dy,  heL2(D)  (17) 

By  |*|,  ||*  ||  we  denote  c(A)  and  L^(A)  norm  respectively. 

Problem  A.  How  to  find  f£eWD  such  that  |f-f£|  <  e? 

Problem  A.  How  to  find  f££WD  such  that  j|f-f£  j|  <  e? 


Let  D  =  meas  D, 


*n(y}  =  [w  ^  expC'2HTN  (t'y,}dtl 


2n+N 


(1  -  ^^-)n  (^y) 

R  TtR 


n  .  M/2 


(is: 


Here  |d|  *  meas  D,  R  >  0  is  such  a  number  that  the  ball  |y|  £  R 
contains  all  vectors  t  -  y,  t,  yeA,  the  origin  in  Rfc  is  placed  in 
the  gravity  center  of  D  so  that 

/D  t  dt  =  0.  (19) 

We  assume  that  A  =  A-  and  D  =  D  are  bounded  simply  connected  and 
connected  domains.  Let 

fn(x)  =  ^A  gn (x-y) f (y) dy  .  (20) 

Theorem  2.  There  exists  n  =  n(e)  such  that  f(x)  is  the  solution 

r\ 

2 

to  problem  A  for  n  >_  n(e)  if  f(x)eC(A)  and  to  problem  B  if  feL  (A). 
Theorem  3.  If  |f(x)|  <  a,  | V f ( x )  ]  <  b,  then 
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If  -  fnl  i 


r(S+i) 

^nT- 


r  (j)  /n 


bR  i  r\  / 1  \ 

—  +0  (— )  as  n 
—  n 


r  (x)  being  Gamma  function. 

Proof  of  theorem  2.  Note  that  g  (y)eW_.  Hence  f  (x)eW  .  Let 


2n+N 


-*■  1  as  n  -*■  », 


(21) 

(22) 


us  prove  that 

[TTT 

/.  (1  -  ly-^1  )n  f(y)dy  -*  f(x)  as  n  -  “ 

ttR  R 

uniformly  in  A  if  f(x)eC(A). 

We  have 

(2n+N)  In  yjjjy  JD  exp  {-  (t,y)}dt  =  (2n+N)  • 

In  [1  2n+N  ^ '  j  D  j  ^ D  tdt^  +  an^  '  an  =  0  (n  ^ 

From  here  and  (19)  we  get  (21).  To  prove  (22)  we  can  assume  with¬ 
out  loss  of  generality  that  R  =  1,  otherwise  we  can  use  scaling. 
For  any  n ,  0<n<0.5we  have 

2,n 


lim 

n-*<»  n<  !  x-y  I  <1 


/  [ 1  -  ix-y |  ]ndy  =  0  , 


lim  (^)N//2  J(  [1  -  |x-y|2]ndy  =  1. 

ix-y|'<n 


(24) 


(25) 


Further 


f  -  f  = 


=  J1  +  J2  H  /a  ^  (y)  "  f(x)]dy  + 


2 ,  n 


+  f(x)  t(2.)N/2  JA  {1  -  I  x-y  |  ^  }ndy  -  1]  , 


2 1  n . 


IJ2I  i  I f ! • I (£)N/2 


U 


(1  -  u2)ndu  -  l|  0  as  n 


<1 


(26) 

(27) 
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|  J,  1  <  (2.)N/2  /  (l-|u|2)n|f  (u+x)-f  (x)  |du  + 

|  U  |  <_£ 

+2|f|(-)N/2  /  ( 1— | u | 2)ndu  (28) 

*  e<|u|<l 

From  (28)  it  follows  that  -►  0  as  n  00 .  We  proved  the  first 

statement  of  theorem  2.  To  prove  the  second  statement  we  must 
prove  that 

2 

||(-Ar) N/2  JA  (1  -  liSzTL.)11  f(y)dy  -  f(x)i|  -  0  as  n  -  ».  (29) 

As  above  we  assume  that  R  =  1.  Let  us  fix  e  >  0  and  find 
4>(x)eC(A)  such  that  ||f-<j>  j|  <  e.  Setting 

Qnf  =  (^)N/2  /A  (1  -  | x-y  S  2 ] n  f(y)dy  (30) 

we  have 

llQnf  -  f II  <  llQn(f  -♦)!!+  -  <M(  +  (If  -  ♦((•  (3D 

Since  |[p|l  <_  C|p|,  C  =  C(A)  we  obtain 

||o  ♦  -  ip  ||  0  as  n  -*■  «.  (32) 

It  remains  to  prove  that 

llQ^II  -,-0asn-*-<=°,<('  =  f-  ,K  (33) 

We  have 

llQ^ll2  =  dx((^)N/2  JA  [1  -  I  x-y  |  2] n  |  ’My)  dy(£)N/2- 

•JA  (1  -  |x-z|2]n|^(z) |dz)  <  i(J)N  J  dudv(l  -  |ui2)n 

A  1  *  |  u  I  <i,|  V i  <1 

{ /A  (|<Mx+u)  |2  +  I  tp  (x+v)  |  2)dx)  <_ 

<  e  (£)N  /  (1  -  |u|2)n(l  -  i v 1 2 ) n  dudv  <  const  e.  (34) 

| U I <1 1  V I <1 

where  const  does  not  depend  on  n. 
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Proof  of  theorem  3.  To  prove  this  theorem  we  must  sharpen  some  of 


the  estimates  given  in  proof  of  theorem  1.  First  we  note  that  an 
in  (23)  can  be  estimated: 


ct  < 

n 1  — 


Jnl  (t,y) |  dt  £ 


2 | D | (2n+N) 

From  (19)  and  (7)  it  follows  that 
1  t  _  r  i  2n+N 


(2n+N) 


(35 


[JD T  eXp  {”  2HTN  (t'y,)dtl 


-  1  =0  (-)  as  n 

n 


(36 


Further 

(~VN/2  Ja  dy  =  1  +  0  (i)  as  n  -  »  (37 

a  R  n 

Now  we  estimate  integral  from  formula  (10)  without  assumption 
R  =  1.  We  have 

J  <  b  (-£*•) N/2  /  (1  -  lii-)n  |  u  |  du  = 

ttR  |u|<R  R 


-  bE  '7>N/2  s» 


C 1  2 .  n  N 

Jn  (1  -  v  )  v 


dv  = 


bRr . 

— m — --  +  o  (  1 

r(|)  /n 


3/2 


)  as  n 


n 


(38 


2tt  7  N 

where  SN  =  f  (N/2)"  is  the  area  the  unit  sphere  in  R"  ,  itegral 

(1  -  v2)nvNdv  =  jB  n  +  1)  ,  B(x,y)  is  beta-function.  We 

can  evaluate  this  integral  as  n  -►  00  either  using  Stirling  formula 
or  Laplace  method  for  asymptotic  evaluation  of  integrals.  As  a 
result  we  get  (20). 
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3.  Approximation  of  derivatives. 

Let  f(x)cC2,  f^eCQ,  Cn  =  Cn[0,l]  is  the  Banach  space  with  the  norm 


|f  ||_  =  I  max  |  f  }  (x)  |  , 


'n 


j=0  0<x<l 


'O' 


|f  -  f6  ||  _<  5. 


Given  f^  how  to  construct  (x)  so  that  e(6) 


llu5  -  f 


as  5-*-0?  This  is  the  problem  of  stable  differentiation.  Let 


’  h_1[f (x+h)  -  f (x)  ]  ,  0  £  x  £  | 


Ahf  = 


h*1[f  (x+|)  -  f(x-j),  »  £  x  <  i  -  | 
h-1[f  (x)  -  f  (x-h) ) ,  1  -  j  £  x  <  1. 


Theorem  4.  Suppose  that  ||f"  ||  £  M,  h(<5)  =  2(6M  ^)^'^»  Then 

ilf’  “  K(5)fll  1  2  ) 1/2  . 

Proof  of  Theorem  4.  We  have  [(A^f.  -  f’[|  £  jjA^  (f  ^  -  f)  ||  + 

+  ijAj^f  -  f’||  £  26h~^  +  0.5Mh.  Minimizing  in  h  >  0  we  find 

h  =  h(6)  =  2  ( 6 M~ ^ ^ ^  and  the  estimate  ||f’  -  K(5)K^  —  2 

Remark  1.  Suppose  that  feC2n_^/  A^f  =  h  ^  £  Aj|nK(x+khn  . 

“  k=-n 

If  coefficients  A^1"^  are  determined  by  the  system 


(39! 


1/2 


k=-n 


l  kPA^n)  =  n«pl#  0  £  p  £  2n-l,  6 


Pi 


_  fo,P*i 


(40) 


then  |  f  —  f  *  (  £  Mnh^n^‘.  We  can  estimate  ||A^n^  f^  -  f '  ||  and  find 
the  optimal  h  =  h(S). 

Suppose  now  that  u  =  f  +  n,  where  u  =  u(t)  is  a  random 
function,  f  is  a  signal  and  n  is  noise.  The  problem  is  to  find  a 
linear  operator  (estimate)  L  such  that  D[Lu  —  f ’ ]  =  min,  D  is  the 


symbol  of  dispersion.  The  solution  of  this  problem  was  obtained 
by  N.  Wiener  for  stationary  random  processes  [23].  This  solution 
is  difficult  to  find  analytically  and  to  apply  in  practice.  We 
give  here  a  quasioptimal  solution  of  this  problem.  The  idea  is  to 
look  for  a  solution  in  a  subset  of  the  set  of  all  linear  estimates. 

We  give  an  estimate  which  is  easy  to  calculate,  easy  to  apply  in 
practice  and  which  has  an  error  of  the  same  order  as  optimal  estimate 
of  f'.  Let  us  look  for  an  estimate  of  the  form 


Lu  =  A,fn)u  =  h-1  7  A*m)  u(t  +  khnf1) 

n  k=-m  K 

where  coefficients  A^n^  are  determined  by  system  (40), 
A^-  +  0.5,  A(J1)  »  0;  A^2)  =  +  1/6, 


i-zV  “  +  4/3 ,  Aq2)  =  0;  A^*  =  +  0.3, 


a|2)  =  +  0.9,  A^=  +  2.25,  Aq3)  =  0;  A^>  =  +  1/70, 

A^}  =  +  16/105,  A^2 5  »  +  4/5,  A^  *  +  16/5,  =  0. 

Parameter  h  is  to  be  found  from  the  condition  D  [Aj^u  -  f '  ] 

Suppose  that  f  and  n  are  uncorrelated,  ,  Rn  are  covariance  func¬ 
tions  of  f  and  n  respectively.  Suppose  that  f  and  n  are  stationary 

2  2 

random  processes  and  n  =  5v,  6  >  0,  D[v]  =  1,  6  =  D[n],  R  =  6  R 

n  v 


We  have 

D[AjJn)  u  -  f  *  ]  =  D[AjJm)  f  -  f  •  1  +  D[Aj[m)n], 


D[Aj!m)h]  =  h_262  l  Afm)A(m)  R  (h(j-p)m'1)  =  X  h~2<S 2 

j ,p=-m  J  P  v 

D[AjJm)  f  -  f’]  <  Mjn  hm~2, 
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where 


«  -  Du2”-1]  (  *  imi 

ra  \k=-m 

Therefore  the  optimal  h  -  h(6)  is  to  be  found  from  the  condition 

_  ,  .  -2,2  .  ,,  .4m-2 
£  =  \n  h  6  +  Um  h  =  min‘ 

Thus  h  =  [A  u"1  (2m-l)-1] (4m) 51/ (2m)c  =  0  (<52"in)  . 
opt  m  m  ' 

These  results  were  used  for  finding  extremal  points  of  random  func¬ 
tions.  If  a  random  signal  u  =  f  +  <5v,  where  v  is  noise,  and 
feC?,  [f”|  <_  M,  is  a  deterministic  function  defined  on  a  segment 
[a,b]  and  unimodal  (it  means  that  f  has  only  one  extremum  t*  in 

the  segment  [a,b]),  then  the  following  algorithm  for  finding  tg 

-1  1/2 

is  convenient  in  practice.  Let  us  take  h  =  2(6M  )  and 

Yj  *  (2h)_1  [ u  ( t j  +  h)  -  u ( t j  -  h)],  j  =  1,  2,...,  t .  ■  a  +  jh. 

If  V4  y.,,  <0  then  we  assume  that  t.  <  t*  <  t..,. 

J 3+1  3  3+1 


.  2m- 1  ,  . 

h  \4m) 


(2m-l)  ! 


4.  Optimal  harmonic  synthesis. 

Suppose  that  the  system  {$n(x)l  forms  an  orthomormal  basis 
2  ® 

of  L  (D)  ,  f  (x)  =  y  a  d>  (x)  ,  and  instead  of  a  the  numbers  c  =  a  +6 

j  n  n  n  n  n  n 

are  known,  where  5  5*"  *  o26  .  Here  the  bar  means  statistical 

n  m  n  nm 

averaging,  6^  is  the  Kronecker  symbol,  star  means  complex  con- 

2  2 

^ugation.  We  assume  that  <  a  ,  n  *  1,2,....  The  problem  is 

to  find  f (x)  with  minimal  error  if  (cn)  are  given. 

We  give  two  approaches  to  the  problem.  The  first  one  con¬ 
sists  in  finding  N  =  N(o)  from  the  condition 
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e  (a)  =  ||f(x)  -  l  C  <j>  (x)  ||2  =  min,  ||*  ||  =  INI  2  (42 

l  n  n  l2(d) 

The  second  approach  consists  in  introducing  some  convergence  multi¬ 
pliers  pn(q)  and  finding  q  =  q(a)  from  the  condition 


ll£  fPn(q)  Cn  "  anJ<^n(x)^2  =  min‘ 


(43! 


We  assume  that 

lanl  <_  An”3,  a  >  0.5.  (44; 

The  answer  to  problem  (42)  is  the  following: 


N(c) 


_l  Va 
(Ao  A)  , e 


mm 


(a) 


< 


2a 

2a-l  ‘ 


(45] 


If  we  take  pR(q)  =  exp  (-qn)  then  the  solution  to  problem  (4  3) 
gives  <30pt  which  is  the  positive  solution  of  equation 

A2  X,  (1+e)  cqc~^  =  2-2exP(~2q)  a  ,  C(x)  5  l  n-x,  ( 

{ 1-exp (-2q)}  n=l 

3+g 

c  -  2  if  a  >  — »  c  =  2a-l-e  if  2a<3+s.  Here  0  <  e  <  2a-l. 

Generalizations,  examples,  applications  to  signal  processing  and 
numerical  results  are  given  in  [33]. 


§8.  Optimal  solution  of  antenna  synthesis  problems. 

1.  Consider  the  problem  of  finding  the  current  distribution  in 
a  linear  antenna  from  the  antenna  pattern.  This  problem  can  be 
reduced  to  integral  equation 

o 

Aj  -  /_£  j  (z)  exp  (ikz)dz  =  f  (k) ,  -  kQ  <_  k  <_  kQ  ,  (1) 
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(1) 


”  vr-;- 


IT  (-£,£} ,  j  (z) 


2 

where  f(k)  is  the  given  radiation  pattern,  j(z)eL  = 
is  the  desired  current  distribution,  j(z),  f(k)  are  scalar  functions 
It  is  quite  clear  that  equation  (1)  has  no  more  than  one  solution 
and  this  equation  has  the  solution  iff  f(k)ew^  =  ^  where 


class  WD  was  defined  in  §7,  section  2. 

f(k)*W£  but  fe(k)eW£,  |Jf  e  -  f  j|  <  e,IHI 


If  the  given  function 


L  (-W 


then  we  consider  j£  corresponding  to  f£  as  an  approximate  solution 
to  synthesis  problem.  As  operator  A  is  compact  small  perturba¬ 
tions  of  f  can  imply  large  variations  of  j .  So  j£  can  change 
much  when  f£  is  changed  a  little.  This  phenomenon  was  discussed 
in  the  literature  (Minkovich-Jakovlev  [1])  in  connection  with  the 
superdirectivity  of  antennas.  From  practical  point  of  view  we 
should  find  a  stable  current  distribution  which  generates  radia¬ 
tion  pattern  close  to  the  desired  pattern  f(k). 

So  we  require  that 

|j|2dz<MQ,  \L_t  ]  j*  (z)  i2dz  <  Mr  (2) 


where  M^,  are  some  constants. 

The  set  of  functions  satisfying  conditions  (2)  we  denote 

2 

by  M^) .  It  is  convex  and  compact  in  L  .  So  we  came  to  the 

problem  of  solving  equation  (1)  under  conditions  (2)  in  the  follow- 

.  2  2 
ing  sense:  given  a  function  g(k)£L  (-k0,kQ)  *  L  we  look  for 

f(k)eW^  such  that 

J  ■  fZao  1 9  00  “  f(k)|2dk  »  min,  g(k)  =  0  for  |k|  >  kQ,  (3) 

f (k)  =  Aj  (see  (1) ) ,  (4 ) 

jeft (Mq/  Mx).  (5) 
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This  optimization  problem  can  be  solved  in  principle  by  non-linear 
programming  methods,  by  direct  methods  of  calculus  of  variations 
and  by  method  of  calculus  of  variations  based  on  the  Euler  equa¬ 
tion.  We  consider  these  possibilities. 

2 .  Let 

g(z)  =  g (k)  exp(-ikz)dk  (6) 

By  Parseval  equality  we  get: 


27  =  /-«  ls<*)  -  j  (2)2dz  =  J  | g  (z) | 2dt  +  |g(z)  -  j (z) | 2  dz 

I  z  j 

=  S(g,£)  +  Jx,  (7) 


where  6  (g,£)  does  not  depend  on  j  ( z ) .  Hence  J  and  reach  their 
minimum  at  the  same  function  and  problem  (3)- (5)  is  equivalent  to 
the  problem: 


{ 


=  min 

jsft (Mq/  Ml) 


(8) 


We  assume  at  first  for  simplicity  that 

j  (l)  =  j  (- 1 ) .  (9) 

In  the  end  we  eliminate  this  assumption.  Let 


_  r  j  exp  (imrz/£)  2  _  v  2  exp  (inirz/.d) 

J  '■z>  L  Jn  >  y  -  A  yn  _  < 

-  “  /2Z  n=-°°  '/TZ 


(10) 


n=-°° 


j  .  ,1  j(2)  exp(-imt2.(_i)  d2 

n  1  /i r 


We  note  that  assumption  (9)  will  be  used  only  once:  it  allows  to 
differentiate  the  Fourier  series  of  j (z)  termwise.  Substituting 
(10)  in  (8)  we  get 
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h=-°° 


3n  -  9, 


l  !:n|2  £  M  0'  X  "2l  jn|2  -  M1  h-  (] 

n=-°°  n=-“  it 

It  is  the  problem  of  convex  programming.  Since  u  (Mg,  M^)  is 

uniformly  convex  and  compact  in  t  and  functional  =  ||g  -  j  || 

is  strictly  convex  problem  (11)  has  a  solution  and  the  solution 

is  unique.  A  numerical  solution  to  problem  (11)  can  be  obtained 

so.  We  fix  N  and  set: 


rlN  =  n  l  l3n  '  V 
n=-N 


\  N  _  N  -  -  »2 

I  I3nl  <  V  „  J  n  ^n1  -M1  7 

n=-N  n=-N  7T 


This  is  a  finite  dimensional  problem  of  convex  programming.  Denote 

(N)  N 

by  j  the  solution  to  problem  (12)  and  show  that  j  ■+  j  as  N  -  00 . 

Here  j  is  the  solution  to  problem  (8)  and  -*■  denotes  convergence 

in  l*.  We  denote  (jN)  »  6N, 


N 

l 

n=-N 


jn  - 

n{  >N 

'  SN  -  dN 

we  get 

lim 

N-»oo 

*  lim  6  , 
N-*-» 

6  <  d. 

lim 

N-*<* 

.N 

D  * 

i  Let  5  *  lim  5N#  5  -  d'  3  *  j 

(  00  ) 

J,  (j  )  ■  5  <_  d.  Hence  by  definition  of  d  we  get  6  *  d  and  by 


uniqueness  of  solution  to  problem  (8)  we  get  j 


j .  Hence 


j  +  j  as  N  ♦  *,  If  geft(Mg,M^)  then  the  solution  to  problem 

(11)  is  trivial:  j  ■  gn  but  it  is  not  interesting  since  in 
practice  gift (Mg, M^).  If  g$ft(Mg,M^)  then  from  geometrical  reasons 
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it  is  clear  that  the  solution  to  problem  (11)  lies  on  the  boundary 
of  the  set  ft(Mg,M^)  and  hence  at  least  in  one  of  inequalities  (11) 
should  be  equality.  If  one  of  inequalities  (11)  is  strict  and  the 
other  is  equality  we  can  omit  the  strict  inequality  and  consider 
the  problem  (11)  with  only  one  restriction  in  the  form  of  equality 
This  relatively  simple  problem  we  do  not  want  to  discuss  and  pass 
over  to  the  more  difficult  case  when  both  inequalities  (11)  are 
equalities.  So  we  get  the  following  problem 


n=-“ 


V  '  •  ~  ,  2 

l  iDn  -  gnl  =  mxn 


2 1  .  ,2 


n=-co 


I  I5nt ‘  -  V  I  . 


n=-° 


M2  =  M1  ~2 
IT 


The  same  is  valid  for  problem  (12),  for  which  we  have 


N 


l  '3n  -  gr 


=  mm 


h*-N 


N 


n*-N 


N 


M 


0'  I  n  ijn 
u  n— N 


2  '  !2  »  M.. 


The  latter  problem  can  be  easily  solved  with  the  help  of  the 
Lagrange  multipliers.  We  set 


$  "  J1N  + 


X  (  l  ljnl2  -  V  +  u(  l  n2^n'2  -  M2} 

n»-N  n*-N 


(13) 


(15) 


rP  *  °'  »0-N<n<_N 


and  write  the  Euler  equations 

3* 

Jn  -n 

System  (16)  can  be  easily  solved: 


1+X+un 


j  ,  -  N  <_  n  £  N. 


(16) 


(17) 
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Parameters  A,  y  should  be  found  from  two  equations 

2 


N 


N 


n2  |gT 


y  - =  Mn ,  l  - Vy  = 

n=-N  (1+A+yn  )  u  n— N  (l+A+yn'V 

Equation  (18)  can  be  easily  solved  by  numerical  methods. 


(18) 


3.  Problem  (13)  can  be  written  so 


J1  =  l3U)  -  g(z)|2dz  =  min 


(19) 


j | 2dz  =  MQ/  | j  ’  (z) i2dz  =  Ml 


Using  calculus  of  variations  we  get  the  Euler  equation  and  the 
natural  boundary  conditions: 

(l+A) j  (z)  -  yj ’ •  (z)  =  g(z) ,  -i  <  z  <  £,  j '  (-£)  »  j  '  (£)  »  0  (20) 

This  equation  can  be  solved  explicitly  and  then  parameters  A,  y 
can  be  found  in  principle  from  equalities  (19).  Problem  (19) 
can  be  also  solved  by  direct  methods  of  calculus  of  variations, 
for  example  by  Ritz  methods.  If  we  set 
N 

I  9 

n=-N 


^  (N) 


exp  (in^rz /£) 

n  STZ 


we  get  system  (16)  for  unknown  coefficients  j  with  restrictions 


(14) 


4.  Here  we  eliminate  assumption  (9).  If  j  (t)  ft  j  (-£)  then  even 
for  a  smooth  in  the  closed  segment  [-£,£]  function  its  Fourier 


coefficients  are  0  (— )  and  so  the  Fourier  coefficients  of  j  * ( z ) 

n  n 


are  not  equal  to  j 


It  is  so  because  the  2£-periodic  continua¬ 


tion  of  the  function  j (z)  is  not  smooth.  To  avoid  this  difficulty 


we  can  use  the  following  orthonormal  basis  in  L  (-1,1): 


*n  *  ^  cos  (£TI  +  n  m  °>  2-*‘ 


(21) 
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Note  that  4>n  are  the  eigenfunctions  of  the  following  problem 

y'  '  +  \y  =  0,  i  z  <  f;  y#  (-O  *  y'  (")  =  0.  (22) 

Fourier  coefficients 

Cn  =  /_  j  (z  )<t>n  (z)dz  (23) 

are  0  (^)  as  n  -►  °°  and  we  can  differentiate  Fourier  series  of 

n  * 

j(z)  according  to  the  system  {«|>  J  toimwise. 

6.  By  the  given  in  this  section  method  one  can  solve  other  antenna 
synthesis  problem,  for  example,  the  problem  of  synthesis  of  the 
spherical  antenna  and  directorial  antenna* 


59.  Stable  solution  of  integral  equation  of  the  first  kind  with 
logarithmic  kernel. 


1.  Consider  the  equation 

Af  =  /  In (r”^ )  f(y)dy  =  g(x),  xeD,  \  =  /„,  (1) 

where  OCR m  is  a  bounded  surface.  Because  Inr  ^  changes  sign,  it 

xy 

is  possible  that  for  some  D  the  homogeneous  equation  (1)  (equation  ( 1q)  ) 

will  have  a  nontrivial  solution.  Let  us  assume  that  (A):  Equation  (l) 

2 

is  solvable  in  H  =  L  (D)  and  equation  (1Q)  has  a  nontrivial  solution  in 
H.  We  give  an  iterative  process  for  calculation  of  solutions  of  equa¬ 
tions  (1),(1q).  The  method. given  holds  for  any  A  semibounded  below 
and  such  that  its  kernel  A(x,y)  >  -k,  k  =  const. 


2.  Let  us  take  a  number  d  >  diam  D.  Then  In  (dr  ^)  >  0  for 

xy 

x,y>D.  Equation  (1)  is  equivalent  to  equation 


3f  ( 2 1  )  ~ 1  l  In(dr^)  f(y)dy  =  h(x)  g(x)  +  lnd  J  f(y)dy  =  g+c(f)  (2) 

Operator  B  is  positive  in  H.  Thus  Bf  =  0  implies  that  f  =  0 .  (3) 


-1, 


Lemma  1.  dim  Ker  A  =  1. 

Proof  of  Lemma  1.  Suppose  f  ,  f^cKer  A. 
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Then  / f^dx  ^0,  j  =  0,  1.  Indeed  if  / fdx  =  0 ,  Af  =  0 ,  then  Bf  =  0 
and  thus  f  =  0.  Let  us  take  b  >  0  such  that  /  (fg  -  bf^)dx  =  0. 
Then  fg  =  bf^.  This  completes  the  proof. 

According  to  Lemma  1  the  problem  Afg  =  0,  /fgdx  =  1  has  a  unique 
solution.  If  equation  (1)  is  solvable  in  M,  it  has  a  solution 
which  satisfies  the  condition 

/ fdx  =0,  (4) 

and  thus  Bf  =  g.  Indeed,  if  Af^  =  g,  then  f  =  f^  -  afg,  where 
a  =  / f^dx,  Af  ^  =  Af  =  Bf  =  g. 

3.  Let  us  look  for  the  solution  of  (1)  which  satisfies  (4).  This 
problem  has  the  unique  solution  for  gsR(A)  and  is  equivalent  to 
the  problem 

Bf  =  g.  (5) 

The  kernel  of  equation  (5)  is  pointwise  positive.  Therefore  the 
theory  given  in  §3  is  applicable  and  we  have  Theorem  1.  Let 

a(x)  =  {/in  (dr'^Jdy}  \  f  =  aij>, 
xy 


*n+l  =  (I  -  V^n  +  g'  *0  =  g' 
where  B.:|>  =  B(aip).  Then  |j^  -  | 


-*■  0  as  n  -*■  ». 


Remark  1.  Iterative  process  (6)  allows  us  to  aclculate  the  solu¬ 
tion  of  (5).  Because  B  >  0  in  H,  there  is  a  possibility  to  use 
iterative  process 

<Pn+1  "  (I  *  kB^n  +  kg,  0  <  k  <  2  | |B  [ |“ 1  (7! 

(see  (4.12)). 


10.  Bibliographical  note 


1.  Much  more  examples  of  ill-posed  problems  can  be  given. 

2.  The  result  given  in  Theorem  1  is  similar  to  some  well-known 
results . 

3.  Theorem  1  was  announced  in  [24] . 

4.  Theorem  1  was  proved  in  [25]. 

5.  The  results  of  this  section  appeared  in  slightly  different 
form  in  [26]  and  were  used  in  [27]— [28] . 

6.  The  theory  outlined  here  was  developed  in  [15] -[19].  The  results 
of  this  section  are  new. 

7.  The  results  of  this  section  were  given  in  [29]— [33] . 

8.  The  results  of  this  section  were  presented  in  [34],  [35]. 

Similar  theory  is  of  interest  in  optics  [29]. 

9.  The  result  of  this  section  is  new. 

The  nonlinear  antenna  synthesis  problem  was  sutided  in  [30]. 

We  do  not  discuss  in  this  paper  some  known  conceptions  in  treating 
of  ill-posed  problems:  Regularization  method  [1],  [3],  [5], 

method  of  quasiinversion  [2],  method  of  quansisolution  [4],  The 
iterative  process  for  solution  of  integral  equation  of  the  first 
kind  with  self adjoint  nonnegative  operator  first  was  given  in  [44] 
(see  also  [5]-[9]).  The  inverse  scattering  problem  was  treated 
in  [11].  In  [14]  the  inverse  problem  of  potential  theory  was 
studied.  In  [37]  the  stability  of  calculations  in  variational 
methods  is  studied.  In  [38],  [39]  inverse  diffraction  problem 

was  studied.  In  [40]  a  method  of  stable  solution  of  some  ill- 
posed  problems  is  given  and  applications  to  partial  differential 
equations  are  considered.  In  [41]  many  inverse  problems,  which 
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are  incorrectly  posed,  are  discussed.  Among  them  there  are  problems 
of  interest  in  geophysics,  optics  and  quantum,  mechanics.  In  [42] 
statistical  approach  to  solution  of  ill-posed  problems  is  discussed. 
Any  linear  solvable  equation  Af  =  g  with  a  bounded  operator  A  on 
a  Hilbert  space  can  be  solved  by  an  iterative  process: 

fn+1  =  (I  -  aA*A)  f  +  ag,  fQeH,  0  <  a  <  2  jjA*A||. 

One  of  the  first  results  in  this  direction  was  given  in  [43].  We 
have  no  possibility  to  discuss  here  the  statistical  approach  to  the 
ill-posed  problems  [42],  [43]. 
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